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Abstract 

We calculate the nucleon-antinucleon static potential in the Skyrme model 
using the product ansatz and including some finite Nc (number of colors) 
corrections. The mid and long range part of the spin-spin and tensor force 
are correctly given in both iso-spin channels while the central interaction has 
insufficient mid-range attraction. This is a well known problem of the product 
ansatz that should be repaired with better Skyrme dynamics. 
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I. INTRODUCTION 



The Skyrme model |l| is an example of what QCD might look like in the classical 
or large number of colors (Nq) limit 0,0]. The dynamics of the SU(2) Skyrme model 
is carried purely by a classical pion field. Hadrons appear as topological solitons in this 
non-linear meson field theory. These are the appropriate degrees of freedom for the non- 
perturbative, long wavelength limit of QCD, and hence for low energy baryon and pion 
physics. The Skyrme model has been applied to the nucleon static properties M and the 
nucleon-nucleon interactions |5],|6],|7]|| wr th reasonable success. In the last few years, nucleon 
annihilation has been investigated from the Skyrme point of view. Sommermann et al. || 
studied the dynamics of ungroomed Skyrmion-AntiSkyrmion (SS) collisions. They found 
that annihilation proceeds quickly with the creation of a coherent pion pulse. This was 



confirmed by Shao, Walet and Amado ||10|| . The notion that annihilation leads to an intense 
coherent pion field burst gives reason for considering annihilation within the classical Skyrme 
approach. The idea, that the annihilation products, pion and other mesons, come from a 
coherent wave of meson fields arising from soliton-antisoliton dynamics, turns out to be 



very fruitful |Tiyi2|| . Experimental data such as annihilation branching ratios among meson 
types and pion charge types from low energy annihilation are well explained with minimal 
parameters \T2\. Furthermore this picture provides a unified view of annihilation in which 
all the channels come from a single process. 

Previous studies of annihilation in the Skyrme context, with the exception of ||, have 
concentrated on the final state mesons. A full account of the process requires a description 
of the initial state nucleon-antinucleon interaction and of the dynamics leading up to anni- 
hilation as well. In this paper, as a first step in that directions, we extend the application 
of the Skyrme model to the interaction of NN in the product ansatz. (We note that the 
energy of SS in the product ansatz was studied for two configurations by Musakhanov and 
Musatov iffi .) Phenomenologically, the NN potential is not as well established as the NN 
potential. At distances less than one fermi, the interaction is dominated by annihilation. 
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However, at larger distances, a meaningful potential can be defined and studied either by 
G-transformation on the NN meson exchange potential or phenomenologically. Here we will 
compare our Skyrme model results to this phenomenology. We will see that at large distance, 
where the product ansatz makes the best sense, the potentials we find agree qualitatively 
and, in most cases, quantitatively with phenomenological interactions. At intermediate and 
short distances, we do less well, but at these distances the product ansatz is not valid. How- 
ever it is still suggestive. To obtain the interaction at intermediate distances, we would need 
to study the full Skyrme dynamics at these distances. For the static SS this is somewhat 
more complex than in the corresponding SS case studied by Walhout and Wambach, ||14|| , 



but is possible and we plan to return to it. The full, time dependent, dynamical SS prob- 
lem is far more difficult than the SS case and is plagued for SS by numerical instabilities 



PJT5[|. For all these reasons, and because this paper is a first step, we begin by exploring the 
interaction in the product ansatz. In order to carry out our study it is necessary to include 
A and A mixing in the potential as was first suggested in the NN case [|T6l,|3||. 



In the context of the product ansatz, we find that the ungroomed SS channel studied 
in is the most attractive channel and that it leads to rapid annihilation. For non-zero 
grooming, we find that the SS interaction can be repulsive. Therefore it seems likely that 
the dynamics in groomed channels could be very different from that exhibited in |§. Since 
the physical nucleon is represented by an average over differently groomed Skyrmions, anni- 
hilation in the nucleon-antinucleon system may proceed more slowly than that of ungroomed 
SS. 

In Sec. 2 we study the interaction energy of SS as a function of separation and relative 
grooming in the product ansatz. We start by studying the very simplified case of Skyrmion 
and groomed antiSkyrmion on top of each other. This is a physically artificial case, but 
it permits analytic evaluation and teaches us something, albeit qualitative, about the de- 
pendence of the SS interaction on grooming. We find that for zero separation and zero 
grooming, the SS system has zero total energy, as we expect. This is the case of complete 
annihilation. However at relative grooming angle of ir, and still zero separation, the product 



ansatz gives a total energy for the SS system of four times the single Skyrmion energy, 
corresponding to a very repulsive interaction. This clearly indicates that the SS interaction 
is a strong function of relative grooming. Next we study, always in the product ansatz, the 
interaction energy of SS at non-zero separation as a function of grooming. We project to the 



nucleon space by the algebraic methods of |17j which also include finite Nc corrections. In 



Sec. 3 we consider the effects of rotational excitations by including intermediate states with 
A and A. We first evaluate the corrections to the NN potential in perturbation theory and 
then study the effect fully by diagonalizing in the space spanned by N, A and corresponding 
antiparticles. Our results are presented in Sec. 4. 

II. THE INTERACTION ENERGY IN THE SS SYSTEM AS A FUNCTION OF 
SEPARATION AND RELATIVE GROOMING 

A. The case of an S on top of a groomed S 

We calculate the energy of the Skyrmion and antiSkyrmion system using the Skyrme 
lagrangian. The density of this lagrangian is given by 

C = fiTr(^CW+) + J^Tr(Q Ml/ Q^+) + £m*Tr(U - 1), (1) 

where U is a unitary SU (2) valued field and 

Qpv = [(9^)11+, [d u U)U + ] . (2) 

The first term in the lagrangian comes from the non-linear cr-model and the second is the 
Skyrme term. The third term is a pion mass term and we take = 139 MeV. We take 



the parameters in the lagrangian to have the values /„. = 93 MeV and e = 4.76 [I4|. These 
values guarantee that the long distance tail of the nucleon-antinucleon interaction will agree 
with phenomenology, by virtue of the Goldberger-Treiman relation. 

We begin by studying the energy in the product ansatz for the case of zero separation. 
We include the first two terms in the lagrangian. The mass term is neglected since it does 



not lead to additional understanding of this simple configuration. Let U = exp(ir • rF(r)) 
be the ungroomed Skyrme SU(2) field. The ungroomed S would be W . The rotation or 
grooming matrix C on S is 

C = cos([3/2)+iT-hsm((3/2). (3) 

This corresponds to a grooming rotation through angle (3 about the n axis. The product 
ansatz with this relative grooming is 

U PA = UCU^Cl (4) 

Note that the energy is a function of only the relative grooming and should be zero with no 
grooming (since for (3 = 0, Up A = 1). 
The energy density is 

S = ~trAA-ltr[A,£i] a (5) 

in Skyrme units (energy in ef n and length in l/(e/„-)). The chiral (left handed) derivative is 

Ci = l4 A diU PA . (6) 

Suppose for the Skyrme chiral angle F(r), the mass contribution to the B = 1 skyrmion 
from the non-linear cr-model term is M 2 and from the Skyrme term M 4 . After some algebra, 
we arrive at the following result for the total energy of the SS product ansatz at relative 
grooming angle (3 and zero separation. 

M S s(/3) = I sin 2 M 2 + f sin 4 (^j M 4 . (7) 

For the profile F(r) which minimizes the B = 1 Skyrme mass, we have M 2 = M 4 = M/2 
from scaling arguments. Here M is the mass of the Skyrmion. We then have 



M s - S {(3)/M = ^sin 2 



l + 2sin 2 | - 



(8) 



which is plotted in Fig. |l|. 



The maximum occurs at (3 = n, the maximum grooming, with the value AM. This 
indicates that the 5* 5" is quite repulsive at this setting. Of course for (3 = or 2tt, the total 
energy is zero. This shows that, even for the artificial case of zero separation, the SS energy 
is very dependent on grooming. 

B. SS at a separation and a relative grooming 

We now study, in the product ansatz, the energy of the SS system as a function of the 
separation and relative grooming. This can only be done numerically. We now use the full 
lagrangian including the finite pion mass term. We put the S and 5* on a 3D lattice with 
the two solitons a distance R apart on the x axis. The product ansatz is 

U PA (r) = U(r-^)CU\r + ^)C\ (9) 

where U (r) is the SU(2) field for a single Skyrmion and C is the grooming matrix. In Skyrme 
units for length (l/ef n = 0.45 fm), the spatial extension of lattice we use is 20 x 10 x 10. 
We evaluate the derivatives of the [/-field by the two-point difference: 

ViU = [U(r + ha) - U(r - ha)] . (10) 

With 64 x 32 x 32 points on the lattice and h = 0.001, we find that, for large separations, 
the total energy is within one percent of twice the single Skyrmion mass. We calculate the 
energy for three interesting configurations: 

1. no grooming (HH) 

2. relative rotation of 7r around x-axis (x-tt) 

3. relative rotation of 7r around z-axis (z-ir) 

The results are shown in Fig. [| Recall that the separation is along the x axis. 

It is instructive to compare the energy of the SS system with the corresponding result 
for the SS system |||6],|7|] . The z-7r grooming is the most attractive configuration for SS, 
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while it is the most repulsive for SS. The x-n grooming, while being the most repulsive 
for SS, is mildly attractive for SS. Finally the ungroomed case is the most attractive for 
SS while it is mildly repulsive for SS. This also leads to speculation about the speed the 
speed of annihilation. In the calculation of ||, the starting configuration is ungroomed and 
annihilation happens very fast. This may reflect the fact that this is the the most attractive 
channel. The physical nucleon is a linear superposition of groomed Skyrmions and perhaps 
in this case annihilation will proceed slower than that seen in ||. 

C. Expansion of the SS energy in the relative grooming variables 

We now turn to an expansion of the energy in the relative grooming variables as a first 
step in obtaining the projection of the SS interaction onto the iViV interaction. We follow 
the methods developed for obtaining the NN interaction from the SS. As in the calculation 
of P,|P7f for SS, the energy for SS can be expanded in the variables c 4 and c • R, with the 
relative grooming matrix C = c^ + ir ■ c and R the vector connecting the centers of the two 
solitons. For the SS, the full expansion is 

V(R, C) = V 1 + V 2 c\ + V 3 (c ■ R) 2 + V 4 ci + V 5 cl(c ■ R) 2 + V 6 (c- R) 4 (11) 

where Vi, i = 1..6 are functions of R. For SS, the symmetry of R — > —R is broken by the 
product ansatz and we need three additional terms for a consistent expansion 

Vgg = V(R, C) + V 7 c 4 (c ■ R) + V 8 c\{c ■ R) + V 9 c 4 {c ■ Rf. (12) 

These terms odd in R are an artifact of the asymmetry of the product ansatz and should 
be discarded. One can use the symmetrized energy (V^g + Vg S )/2 to extract V\ to Vq, since 
the Vj to Vg terms drop out in this combination. 

The six terms in (11) can be expressed in terms of operators in the baryon space using 
the algebraic methods introduced in ||17|| . One quantizes each Skyrmion with a w(4) algebra 
and then the relevant operators and baryon states are easily constructed in terms of the 
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operators of those algebras. The method was developed in [[H] for the NN system, but 
since each Skyrmion gets its own algebra, the method can be taken over without alteration 
to the NN system. The SS or SS interaction can be expanded in terms of three operators, 
the identity and the operators W and Z given by 



W = T^/N 2 C , 



pi P3 



(13) 



Here a and (3 label the two different set of bosons, used to realize the w(4) algebras. T is 
a one-body operator with spin 1 and isospin 1. The semiclassical (large- Nc) limit of these 



operators can be given in terms of R and C as [17 



W cl = 3c 2 - c\ 
Z cl = 6c- R- 2c 2 . 



(14) 



We can therefore expand the interaction in W and Z, as an alternative to Eq. (11). 



V(R, C) = Vl {R) + v 2 (R)W cl + v 3 (R)Z cl + v 4 (R)W* + v 5 (R)W cl Z cl + v 6 (R)Z* (15) 



The relations between Vi and Vi can be found in Eq. (24) of JT^]. The advantage of the 
algebraic method is it allows us to study both the large Nc limit, and to include finite Nc 
effects explicitly in a systematic way. It also makes taking baryon matrix elements quite 
easy. As in the SS case, we find for SS that the terms quadratic in Z and W are quite 
small, and so we neglect them. Hence we can write 



V = v 1 + v 2 W + v 3 Z. 



(16) 



We can use the algebraic methods of |17j to take the NN matrix element of our interac- 
tion. Keeping only the leading terms we find 



V = V c + K(<Ti • <r 2 )(n • r 2 ) + V t [Z(tr x ■ R)(tr 3 ■ R) - a ± ■ tr 3 ](n • r 2 ) 



(17) 



with 



V c = v u V. = *pL, V t = V -^L. (18) 

Here Pn is a finite Nc correction factor, Pjy = 1 + 2/Nc- This gives the nucleons only 
projection of the SS interaction. To obtain the full phenomenological interaction it is 
necessary to include the effects of A and A admixtures that become important as the 
baryons approach each other. We now turn to those admixtures. 

III. ADIABATIC INTERACTION 

The NN potential in Eq. ( |17|) is calculated by projecting Eq. (|16D to the nucleon degrees 
of freedom only. This is certainly the correct procedure for large separation. However as 
the Skyrmion and antiSkyrmion approach, they can deform. In terms of the baryon degrees 
of freedom and Nc = 3 that means excitation of the A and A intermediate states. All 
that is required to define a NN interaction is that the particles be NN asymptotically. 
They may deform or excite as they wish as they interact. In the NN case we saw that 
this intermediate excitation plays a significant role in the intermediate range attraction 
||. For neutral atoms, a corresponding virtual excitation process leads to the attractive 
Van der Waals force at large distance. For Skyrmions, beside this state mixing, there is 
a dynamical distortion that goes beyond the product ansatz. This is a crucial part of the 



NN interaction |14 |, but the corresponding SS distortion is beyond the scope of this paper. 
For the NN system this distortion, coupled with the state mixing, is crucial for getting the 
mid-range attraction. For the nucleon- antinucleon system we expect similar enhancements 
of attraction coming from the distortion. The effects of distortion and state mixing both 
come in at distances where the product ansatz can be expected to fail. Thus our short and 
mid-range results with the product ansatz, even including state mixing, should be taken as 
only indicative and not as the final word. 

As in the NN system, we first study the effect of mixing A and A on the energy 
perturbatively and then use the Born-Oppenheimer method to consider the effect exactly in 
the limited subspace. 
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A. Perturbation Theory 



We first include the effects of the intermediate states, iVA, AN and A A on the nucleon 
antinucleon potential perturbatively. Since we are using separate it(4) algebras for each 
Skyrmion, the results for iViV in || can be carried over to the NN problem, with Eq. (15) 
in II being the perturbation correction. 



Vprp 



Q 



s 



(v 2 2 + 2vi) 



+ (3cr 1 • iter 2 ■ R - a 1 ■ er 2 



.16 



(19) 



18 Q^o T + ( 81 ^v 162 ^) P i T 



V2V3) \ ■ 



Here 6 is the N-A mass difference, Pf is a projection operator onto isospin T, and Qn 
is another finite Nc correction factor with the value Qn = J (1 — 1/AT C )(1 + 5/Nc)- This 



expresses the leading order correction from state mixing to the iViV interaction of Eq. ([Ti 
in terms of the SS terms of Eq. (|i6|). Recall that unlike the work in || we are here using 
the product ansatz to calculate the SS interaction rather than a full dynamical scheme. 



B. Diagonalization 

We now turn to a full diagonalization of the interaction in the NN, NA, AN and A A 
space. This is the Born-Oppenheimer approximation, and it is valid for Nq = 3. There are 
three energy scales or time scales in the problem. The fastest or highest energy scale comes 
in rearrangements of the pion field itself. These we are modeling using the product ansatz 
and correspond to energies on the scale of baryon masses. The intermediate scale is set by 
the N A energy difference. This is an order 1/Nq effect. Finally the iViV interaction is the 
smallest energy scale and it is determined by the matrix diagonalization. 

We first need the matrix element of the potential Eq. ( pp in the space of NN, NA, AN 
and AA in the angular momentum coupled form. This has been calculated in Eq. (22) of 
P| for the baryon-baryon case and the formula remains valid for baryon-antibaryons. 
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(hI 2 LSJT\ v \I[I' 2 L'S'JT) 

= v^ss'Sll'Si^S^ 
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The relevant reduced matrix elements are 



(N\\T^\\N) = -10, 
(A||T( n )||A) = -20, 
(iV||T (11) ||A) = -8^. 

The matrix element of the kinetic part are taken to be very simple, 



(21) 
(22) 
(23) 



{hI 2 LSJT\K\I[I' 2 L'S'JT) = 8 hI[ 5 l2l , 2 5 LL ,5 S s> (s [h + + ^ 



l(l + iy 

h R\ 



(24) 



Here M hh is the reduced mass, M h M h / ( M h +M h ) , with M 1/2 = 932 MeV and M 3/2 = 1232 
MeV. The mass difference is 5M = M 3/2 - M 1/2 . 

For the purpose of comparison, we parametrize the full NN interaction by 



V nn = V c + Vs ^■<r I + V/ alojQlURj - 5a). 



(25) 



The potentials have explicit isospin dependence due to the mixing with states of A. To 
determine the adiabatic potential for NN, we start at large R where we have nucleons only. 
As we move to smaller distance, we diagonalize the K + V matrix and follow continuously 
the eigenvalue corresponding to the NN channel. We then subtract the expectation value 
of K to obtain the adiabatic interaction energy as a function of R. 
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We first consider the case T = 0. For J 71 = 0~ (note that nucleon and antinucleon have 
opposite intrinsic parity), we have three channels \NNL = 5 = 0), |AAL = 5 = 0) and 
|AAL = 25 = 2). The lowest eigenvalue of the Hamiltonian K + V should be identified 
with 

(L — S — 0\V%=°\L = 5 = 0) = V? - 3V S °. (26) 

For J w = + , there are three channels \NNL = 15 = 1), |AAL = 1 5 = 1) and |AAL = 
3 5 = 3). The lowest eigenvalue should be equated to 

(L = 5 = 0\Vffi\L = S = 0) = K° + K° - 4Vf (27) 

We consider one more set of states with J 71 " = 1~ and there are six channels: |iViVL = 
5=1), |iViVL = 25=1), |AA01), |AA21), |AA23) and |AA43). The matrix element 
to identify the lowest eigenvalue with is 

(L = 5 = 1\V^ °\L = 5 = 1) = V c ° + V s °. (28) 

From these three linear combinations of V c , V s and V t in Eq. ( ^6|) to Eq. (^8|) , the potentials 
in Eq. (p5|) are easily solved for T = 0. A similar calculation applies for T — 1, except that 
now iVA and AA^ channels also appear. 



IV. RESULT AND DISCUSSION 

For each total isospin T = 0, 1, we calculate V^. T , V^ T and V^ T in Eq. (p5|) as outlined 
in the last section. We then compare these result with the phenomenological potentials of 
Brian-Phillips |18| and of the Nijmegen 0] group. These are potentials based on meson 



exchange at large distances and phenomenology, including an absorptive part to model 
annihilation, at small distances. The meson exchange part of these potentials for NN is 
obtained from the corresponding NN potentials by G-parity transform - the contribution 
of a particular meson for NN is equal to its part in Vnn multiplied by the meson's G- 
parity. We only compare with the scalar, tensor and spin-spin parts of the potentials. The 
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spin-orbit force is of higher order in 1/Nc and we have not calculated it in the Skyrme 
picture. We should note that various cutoffs are used in the Brian-Phillips, Nijmegen and 
other similar potentials. As a result at distance 1 fm or less the strength of the potentials 
can be significantly different from their meson exchange value. In addition, at distance 
less than 1 fm, the interaction is dominated by the absorptive potential of order 1 GeV. 
Furthermore, at these short distances the entire static Skyrme approach, to say nothing 
of the product ansatz, is no longer meaningful. Hence we should not place any faith on 
comparisons of our results with the phenomenological potentials around 1 fm or less. At 
intermediate distances, between 1 fm and 2 fm, the results from our Skyrme approach to 
the AW interaction suggest that the product ansatz is on the right track, but that careful 
comparison with phenomenological potentials requires a more complete calculation of the 
Skyrme dynamics. In particular we expect the product ansatz to underestimate mid-range 
attraction, as it does for AW. This is basically a consequence of the variational theorem. 
With these thoughts in mind, let us turn to our results. 

Figure 3 shows our results for the T = part of the central potential. We see that the 
effects of A mixing either in perturbation theory or full Born-Oppenheimer diagonalization 
is significant, but still does not begin to agree with the strong central, mid-range attraction 
seen in the phenomenological potentials. This is the fault of the product ansatz we referred 
to above. It will be important to see if complete Skyrme calculations can repair this fault. 
Figure 4 shows the T = 1 central potential. Here the effects of A mixing are more striking 
since for T = 1 single A intermediate states are permitted. Now we do find some central 
attraction, but not as much as is seen phenomenologically. Note that where they differ, the 
full diagonalization result has superior credentials to the perturbation theory result, and it 
is the diagonalization result that is too weak. Again we must await full Skyrme calculation 
in this channel. Figure 5 shows the T = spin-spin part of the interaction. Except at 
the smallest distances, the results are very satisfactory. Note that this is remarkable, since 
the spin-spin interaction is very weak (note the scale in Figure 5), and hence arises from 
cancellation of much larger terms. We believe it is significant that the Skyrme picture 
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can reproduce this scale and even the correct sign. The effects of cancellations are even 
more striking in Figure 6 which shows the T = 1 spin-spin interaction. Again the order of 
magnitudes are correct while the sign depends sensitively on how we do the A mixing. Note 
that the phenomenological potentials are consistent with zero. Finally in Figures 7 and 8 
we show the T = and T = 1 tensor potentials. All calculations of these agree since they 
are dominated by one-pion exchange. (Recall that the Skyrme picture in the product ansatz 
gets one pion exchange right.) Hence except for the central attraction, the product ansatz 
gives a credible account of the nucleon-antinucleon potential, and we understand how the 
product ansatz fails for the central attraction. Note that there are no free parameters on 
our calculation. 

We have shown that the Skyrme picture with the product ansatz is a reasonable first 
step to obtaining the real part of the nucleon-antinucleon interaction. We also understand 
how doing the Skyrme dynamics better can repair the lack of central attraction we find 
here. Hence the next step is to do that dynamics. Then combing this picture of nucleon- 
antinucleon interactions in the entrance channel based on the Skyrme model with our previ- 
ous work on annihilation channels described by this model, we hope to have a unified picture 
of annihilation based on the large N c , QCD inspired Skyrme picture. 
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FIGURES 
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FIG. 1. Mass of the SS at zero separation and with relative grooming angle (3. Ms is the 
Skyrmion mass. 
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FIG. 2. Total energy of the S-antiS system as a function of separation for the configuration HH 
(dashed line), x-ir (dash-dotted) and z-ir (solid) in units of the Skyrme mass. Note the horizontal 
line is twice the Skyrmion mass. The maximum value of the energy at zero-separation is four times 
the Skyrmion mass, as we derived in the analytical result. 




0.5 1.0 1.5 2.0 

R (fm) 

FIG. 3. Central potential as a function of R in the region 0.5-2 fm for the T = chan- 
nels. The solid line gives the nucleons only result from the product ansatz. The short dashed 
line is the result of the state mixing in perturbation theory and the long-dashed line of the full 
Born-Oppenheimer diagonalization. The meson exchange potentials are shown by the dash-dotted 
line for Bryan-Phillips potential |[^] and by the dotted line for the Nijmegen potential 
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FIG. 4. Central potential, same as in Fig. 3 but for T = 1. 
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FIG. 5. The spin dependent potential V s as a function of R in the region 1-3 fm for T = 0. 
Labeling of curves is the same as in Fig. 3. 
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FIG. 6. Spin-dependent potential, same as Fig. 5 but for T = 1. 
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FIG. 8. Tensor potential, same as in Fig. 7 but for T = 1. 
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